Bounds on long-lived charged massive particles from Big Bang nucleosynthesis 
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The Big Bang nucleosynthesis (BBN) in the presence of charged massive particles (CHAMPs) is 
studied in detail. All currently known effects due to the existence of bound states between CHAMPs 
and nuclei, including possible late-time destruction of ^Li and ^Li are included. The study sets 
conservative bounds on CHAMP abundances in the decay time range 3 x 10^ sec ^; Tx ^ 10^'^sec. ft 
is stressed that the production of ®Li at early times T ~ lOkeV is overestimated by a factor ~ 10 
when the approximation of the Saha equation for the *He bound state fraction is utilised. To 
obtain conservative limits on the abundance of CHAMPs, a Monte-Carlo analysis with ~ 3 x 10® 
independent BBN runs, varying reaction rates of nineteen different reactions, is performed (see 
attached erratum, however). The analysis yields the surprising result that except for small areas 
in the particle parameter space conservative constraints on the abundance of decaying charged 
particles are currently very close to those of neutral particles. It is shown that, in case the rates of a 
number of heretofore unconsidered reactions may be determined reliably in future, it is conceivable 
that the limit on CHAMPs in the early Universe could be tightened by orders of magnitude. An 
ERRATUM gives limits on primordial CHAMP densities when the by Ref. [2^ recently more 
accurately determined CHAMP reaction rates are employed. 
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I. INTRODUCTION 

Big Bang nucleosynthesis (BBN) has proven itself as a 
powerful tool in constraining the conditions of the early 
Universe and physics beyond the standard model. Thus 
bounds on a variety of hypothesis have been derived in- 
cluding, for example, modifications of gravity, baryon in- 
homogeneity, matter-antimatter domains, non-zero lep- 
ton chemical potentials, and relic decaying particles. It 
has been recently realized that BBN may also place 
bounds on the abundance of charged, weak-scale mass 
particles existing during and after BBN. Though it had 
been noted already earlier that negatively charged, weak- 
scale mass particles (CHAMPs) form bound-states with 
positively charged nuclei towards the end of BBN , it 
has only recently been put forward that this may impact 
considerably the light-element yields synthesized during 
BBN dii. 

Here the most important proposed change is due to a 
catalysm of reactions such as ^H + '^He — > ^Li -I-7 0. 
Other less important modifications, concerning con- 
straints on CHAMPs, had also been considered Q. 
Being of quadrupole (E2) nature the ^H -I- ''He reaction 
has a very small rate (S-factor: 10~^MeV barn) thus 
yielding typically very little ^Li/-'H'^ 10^^''' in standard 
BBN. When the helium-nuclei is in a bound state (^He- 
X^) the above reaction may be replaced by its photonless 
analogue: ^H -I- (^He-X~) — > ^Li +X~, with a cross sec- 
tion estimated orders of magnitude larger than that for 
the standard BBN ^Li-synthesizing process. Initial esti- 
mates for this enhancemant factor were given at around 
6 X 10'' 0, 0]- It was argued that since ~ 1 of all are 
within bound states with '*He at temperatures T^8 keV 
(cf. Fig. [3]) very small abundances of CHAMPs present 
at t « 10^ sec in the early Universe could already over- 
produce the ^Li isotope with respect to observations Q. 
Assuming a (too restrictive) ^Li/^H ^ 2 x 10~^^ con- 



straint, bounds as strong as rix- /s< 2.5 X 10^1^ the 
CHAMP-to-entropy ratio, were derived. These bounds 
were subsequently weakened by one order of magnitude 
when a more proper evaluation of the rate for the ^H 
+ (*He-X~) ^Li +X~ process was derived. Such 
bounds have now been utilised by a number of groups to 
constrain, for example, abundances of supersymmetric 
staus f s in the early Universe 0, 0, H, Hi El • 

Recently, I have shown that there are several changes 
to the "naive" picture of synthesis of ®Li in the presence 
of CHAMPs [l|]. First, ("He-X") bound states may be 
destroyed during the electromagnetic cascades induced 
by the decay of CHAMPs, rendering ^Li production in 
some parts of parameter space much less efficient. More 
importantly, when CHAMPs are still present at times 
t ^ lO^sec Big Bang nucleosynthesis enters a second phase 
of Coulomb-unsupressed reactions on bound states be- 
tween charge Z = 1 nuclei and a CHAMP. It has been 
shown that reactions such as ^Li -I- (^H-X^)^ *Hc + 
^He +X-, '^Be + (iH-AT")^ ^B +X-, etc, are capable 
of completely destroying any priorly synthesized ^Li and 
^Li. This is possible in particular at somewhat higher 
CHAMP-to-entropy Yx = nx/s ratios. In order to esti- 
mate the efficiency of such destruction, not only was it re- 
quired to estimate the cross sections for such ^Li and ^Li 
destroying reactions, but also those of CHAMP exchange 
reactions such a.s {^ll-X~)+'^E^'^E + (^B.-X-), capable 
of significantly reducing the ^H bound state fraction. Al- 
together nineteen reactions of significant importance for 
the late-time nucleosynthesis t ^ lO^sec have been identi- 
fied. The rates for all these reactions were determined in 
the Born approximation. Concerning details on the BBN 
with CHAMPs at late times, the importance of particu- 
lar reactions, and their evaluation, the reader is referred 
to the original paper [l^]- Unfortunately, the Born ap- 
proximation is not a particular good approximation for 
determining rates of these nineteen important reactions 
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FIG. 1: Limits on the primordial CHAMP-to-entropy ratio 
= nx/s (with nx-/s = Y^^/'i) for CHAMPs with inter- 
mediate life times. Shown are constraint lines for CHAMPs 
of mass Mx = 1 TeV and a variety of hadronic branching ra- 
tios Bh = 10~^ — 1, as labeled in the figure. Solid (red) lines 
correspond to the conservativelimit ^Li/^Li < 0.66, whereas 
dashed (blue) lines correspond to ®Li/^Li < 0.1. It is seen 
that only for CHAMPs with Bu ^ 10"^ the effects of bound 
states become important. For smaller decay times Tx the lim- 
its on CHAMP abundances are virtually identical to those on 
the abundance of neutral relic decaying particles p^ . 



for CHAMP BBN, leaving significant uncertainty in the 
BBN yields with late decaying > lO'^sec CHAMPs. 

Given the above, it seems very premature to rule out 
CHAMPs in the early Universe simply by their produc- 
tion of ^Li at T « 8keV. In this letter constraints on 
the abundances on CHAMPs are derived which take full 
account of all the above mentioned extra physics priorly 
neglected. Here constraints will be placed for two differ- 
ent decay time regimes, for 3 x lO^sec ^Tx^ 5 x lO^sec in 
Section 2, where the important rates are relativily well 
known, and for 5 x lO^sec ^ ^ lO^^sec in Section 3, 
where a Monte-Carlo analysis is employed to derive con- 
servative limits. It will be seen that constraints change 
by large factors with respect to those priorly given, par- 
ticularly for long CHAMP decay times. 



FIG. 2: Same as Fig. [T] but for = 100 GeV. 
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FIG. 3: Fraction of CHAMPs fx which are bound to *He as 
a function of temperature for (a) the approximation by the 
Saha equation, thin-dashed (green) and (b) full numerical in- 
tegration of the rate equation, thin-solid (red). Also shown 
is the product of fx with the rate (av) for the ®Li producing 
reaction ^H(*He— X~ , X~)^Li, simply denoted as "rate" and 
in arbritrary units, for both cases (a) thick-dashed (green) 
and (b) thick-solid (red). The figure illustrates that ®Li pro- 
duction at r ~ 8 keV due to bound states is overestimated by 
a factor ~ 10 when the approximation of the Saha equation 
is utilised. 



II. CONSTRAINTS ON CHAMPS WITH 
INTERMEDIATE LIFE TIMES 

In the next two sections constraints from BBN on 
the existence of CHAMPs in the early Universe are pre- 
sented. Such constraints get increasingly more uncer- 
tain as the life time of a CHAMP increases. For 
Ta; $ 3 X lO^sec CHAMPs have almost no impact on BBN 
beyond those of their injection of electromagnetically and 
hadronically interacting particles during their decay [is'] . 
Such constraints have already been discussed in detail in 
the literature and the reader is referred to, for example, 
Ref. for details. For 3 x lO^sec ^ r^; < 5 x lO^sec con- 



straints are still fairly reliable and depend mostly on the 
^H('*He— AT^, Ar^)^Li rate. Since this one has been de- 
termined beyond the Born approximation and is most 
likely known to within a factor three, limits should also 
be known up to such a factor. It has been found numer- 
ically that other reactions, such as'^H(^He— X~)^Li 
and ■^He(^He— X^, Ar^)''Be play less of a role in setting 
constraints at early times, even in case their rates signif- 
icantly exceed those determined in the Born approxima- 
tion. 

When deriving constraints the following conservative 
observationally determined limits on the light element 
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abundances are adopted: 

Yp < 0.258 

1.2 X 10^'^ < ^H/H < 5.3 X 10"'^ 

^He/^H < 1.52 (1) 
8.5 X 10"" < ^Li/H < 5 X 10^1° 
^Li/^Li < 0.66 (0.1) 

Here Yp denotes the helium mass fraction. The observa- 
tions behind these hmits are discussed in further detail 
in Ref. [l3|. It should be noted here that the frequently 
used bound ^Li/^H ^ 2 x 10~^^ is too stringent. ^Li 
may be destroyed during the life time of a Population 11 
star. In fact, if the current discrepancy between stan- 
dard BBN predicted ^Li/^H « 4 - 5 x lO^^" and Pop 
II star observed ^Li/^H w 1 — 2.5 x 10~^° is resolved by 
factor 2 — 3 stellar ^Li destruction, as claimed for ex- 
ample in Ref. il5], than ^Li is destroyed by at least the 
same factor. This would yield an upper limit close to 
^Li/iH <4 X 10"", or ^Uf^Li <0.1. However, since 
^Li is more fragile than ^Li it may, in principle, be de- 
stroyed by much larger factors than ^Li, as shown in a 
number of stellar evolution studies [l^. A conservative 
limit of ^Li/^Li ^0.66 (corresponding approximately to 
^Li < 2.7 X 10"^°) was therefore apphed in Ref. 14]. 

In Fig. [Hand Fig. [2] constraints on decaying CHAMPs 
with intermediate life times are presented. In both fig- 
ures results are shown for a variety of hadronic branching 
ratios Bh, with Fig. [T] showing results for Mx = 1 TeV, 
and Fig. [5] for ~ 100 GeV. In order to derive these 
constraints the rate for ^H^He— Ar")®Li as given 
in Ref. Q has been utilised. It is seen that for large 
hadronic branching ratio {Bh ^ 0.01 for = 1 TeV and 
5/1^0.1 for Mx = 100 GeV) constraints depend almost 
linearly on B^, and are not different from those for neu- 
tral particles. This seems somewhat surprising, due to 
the advocated power 0, of helium-CHAMP bound 
states to produce ^Li. Nevertheless, it is known that ®Li 
is also produced abundantly by hadronic decays during 
that time [ij, [l3|, and this ^Li source is more impor- 
tant than that of '^li{^lle-X~' , X'fU at large Bh. The 
efficiency of ^H(*He— AT", X~)^Li has in any case been 
overestimated in some papers [1, (T^). ^Li production 
here is given by a convolution of two exponentials, (a) 
the exponentially rising ^He— X~ bound state fraction 
at T lOkeV, and (b) the due to a Coulomb barrier 
exponentially decreasing reaction rate. This leads to the 
bulk of the ^Li production in a very narrow temperature 
interval 8keV>T>6keV. Refs. ^ and fTo"! assumed 
the applicability of the Saha equation for the "'He bound 
state fraction (though earlier studies Refs. [3, d, 0, H 13 
did not). Numerical integration shows that the forma- 
tion of appreciable bound state fractions is slightly de- 
layed when compared to the Saha equation, since the 
recombination rate is of the same order as the Hubble 
rate. This may be seen in Fig. [3] This slight difference 
results than in approximately one order of magnitude 
less ^Li production, due to the convolution of exponen- 
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FIG. 4: Liklihood areas in the CHAMP-to-entropy ratio Yx - 
CHAMP life time tx (in seconds) plane for bound-state BBN 
to obey the observational constraints on light element abun- 
dances Eqs. 1. Results of a Monte-Carlo analysis varying 
nineteen ill-determined reaction rates randomly (see text for 
details) which significantly impact BBN yields at late times 
Tx ^ lO^sec are shown. From bottom to top, the areas show 
likhhoods: > 99% lightest shade (yellow), 95% -99% (green), 
80% -95% (purple), 20% -80% (red), 5% -20% (light-blue), 
1% — 5% (black), and < 1% (dark-blue), respectively, for 
CHAMP BBN to respect observational constraints. No ef- 
fects of electromagnetic- and hadronic- cascades due to the 
CHAMP decay have been taken into account (cf. Fig. 
The Monte- Carlo analysis presented in this figure employs 
~ 1.5 X lO'' independent BBN calculations. 

tials. Very recently, the same observation has also been 
made in Ref. jlSj utilising the set of Boltzman equations 
relevant for ^Li pro duction via ''He— AT" bound states 
as given in Ref. [l^. Conservative limits on CHAMPs 
in the intermediate decay time interval are thus weaker 
than initially thought. Here the weakening is due to a 
lower rate |9| , the failure of the Saha equation, and a too 
restrictive upper limit on the ''Li abundance. It is noted 
here that the photodisintegration of (^He— X~) bound 
states as noted in Ref. ^12i] is comparatively unimportant 
at small Y^, such that further weakening of the limit on 
CHAMPs does not result. In any case, when Bh ^ 10"^ 
bound state induced production of ^Li becomes dominant 
over hadronic production. Since this is the case for, for 
example, supersymmetric staus, bound state effects thus 
still remain very constraining in particular scenarios ^] . 

III. CONSTRAINTS ON CHAMPS WITH LONG 
LIFE TIMES 

When the hfe times exceeds ^ 5 x lO^sec it be- 
comes substantially more difficult to place reliable limits. 
This is due to a large number of CHAMP-induced reac- 
tions becoming important at T ^ 1 keV, in particular all 
Coulomb-unsupressed nuclear reactions shown in Table 
11 and Table HI of Ref. [H], as well as the CHAMP ex- 
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FIG. 5: As Fig.|4]but including the effects of electromagnetic 
cascades during X-decay, assuming that a fraction /em ~ 1 
of the particles rest mass (taken mx ~ 100 GeV) is converted 
into electromagnetically interacting particles. The hadronic 
branching ratio was set to Bh = 0. The dashed line shows 
the analogous limit for neutral relics. 




FIG. 6: Isocontours of ''Li/'^Li = 0.66 solid (red), and ''Li/'^Li 
= 0.1 dashed (blue), in the CHAMP-to-entropy - CHAMP 
life time plane for rates determined in the unreliable Born 
approximation p^ . The present figure is not supposed to 
be utilised for limiting CHAMP abundances. It rather is for 
illustrative purpose, showing the possible importance of the 
*He + (^H—X~) — > X~ +^Li reaction to further tighten limits 
on long-lived CHAMPS by orders of magnitude due to late- 
time ®Li production. A hadronic branching ratio of Bh ~ 
has been assumed. 



change reactions shown in Table IV of that paper. This 
comprises a total number of nineteen reactions. Though 
all rates have been determined numerically in the Born 
approximation in Ref . [l^ , as the Born approximation is 
likely to fail badly, results become uncertain. In order to 
still arrive at a reliable result one is thus forced to per- 
form a Monte-Carlo analysis, varying all ill-determined 
reaction rates within conservative ranges. This has been 
done in the present paper. In particular, the Born ap- 



proximation values of the rates given in Ref. (shown 
in Figs. 3 and 4, as well as in Table III and IV of that pa- 
per), have been taken as benchmarks. For each reaction 
a random generator determined a factor fi with which 
the benchmark rate was multiplied. These factors where 
generated with a probability distribution flat in logarith- 
mic space, and between values l//f"* < /i < /f"*- For 
the reaction-rate dependent conservatively chosen /™* 
the reader is referred to Table VII of Ref. ■ For each 
point in parameter space, i.e. for and t^, this proce- 
dure was repeated a 1000 times in order to arrive with 
one thousand different randomly chosen sets for the 19 
ill-determined reaction rates. For each realization of reac- 
tion rates an indpendent BBN calculation was then per- 
formed and compared to the observational constraints. 

One may wonder if 1000 realizations for the reaction 
rates are actually sufficient for sampling the, a priori, 
complicated probability space. After all even only adopt- 
ing two values for each reaction rate, a large rate and a 
small rate, already yields 2^^ « 5 x 10^ different pos- 
sibilities for sets of reaction rates. For a few individual 
points in and t^, I have therefore generated 10^ differ- 
ent realizations of rate combinations. Comparison with 
the results by only a 1000 realizations shows that the 
simulation with a 1000 realizations may be trusted ap- 
proximately to the 1% level (i.e. < 10 BBN runs passing 
or failing observational constraints). This is true since at 
an individual point in the Y^-Tx plane, results mostly only 
depend on a number '--^ 4 — 6 of rates, with all other rates 
being less important. Which rates are most important for 
the BBN yields, however, depends on the location in the 
Yx-Tx plane. For example, at large Y^ (and t^) results de- 
pend sensitively on the reactions: '^II('*IIe— X^yLi, 
6Li(iH-X-,X-)4Hc + ^He, and {^}i-X-)+ ^R^^B + 
(^R-X-) (and to a lesser degree "^Bci^R-X- , X-)^B, 
^Ue{^}i-X-,X-yU, and {^ll-X-)+ ^He^^H + (^He- 
X^ j), whereas for small Y^ (still large Tx) they depend 
mostly on reactions: ^}ie(^R-X- , X-fLi, (^H-A:-)-!- 
2H^1H + (^U-X-) and {^B-X-)+ ^Be^^U + (^He- 
X-) (and to a lesser degree {^B-X-)+ ^Bc^^B + (^He- 
X-)). 

The results of the present Monte-Carlo analysis are 
shown in Figs. [4] and \5\ Shown are the liklihood ar- 
eas in the Yx-tx plane by different shading (color- 
ing) that > 99%, 95% - 99%,80% - 95%,2G% - 80%, 
5% - 20%, 1% - 5%, and < 1% of all randomly generated 
models at the same Yx and tx obey the observational 
constraints. Note that the CHAMP-to-entropy ratio Yx 
is easily converted to Qx, the fractional contribution of 
the AT-particle to the present critical density, if it had 
not decayed, via Qxh'^ = 2.73 x 10" (A^/l TeV), 
where h is the present Hubble constant in units 100 
kms^^Mpc^^. Note also, that in all constraint figures 
Yx denotes the total CHAMP-to-entropy ratio, assuming 
that only half of all CHAMPs are negatively charged. 
Whereas Fig. 2] shows results when the CHAMP decay is 
not associated with electromagnetic- or hadronic- energy 
release. Fig. [S] assumes a 100% electromagnetic CHAMP 
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decay and the associated electromagnetic cascade nucle- 
osynthesis. It is evident from Fig. [1] that, when only 
effects of bound states are taken into account but not 
injection of energy during the decay, the probability dis- 
tribution is extremely flat. In particular, whereas one 
still finds a small fraction of models which fail for Yx 
as low as ~ 10~^^ only at Yx ^ 10~^^ one may exclude 
CHAMP BBN at the > 99% confidence level. Conser- 
vatively, and in the absence of more reliable rates, even 
Yx as large as 10~^^ may thus not reliably be ruled out. 
This is five orders of magnitude less stringent than the 
initial claim Note, however, that such large Yx may 
only be (conservatively) acceptable for invisible, mass- 
degenerate, or stable CHAMPs, associated with very lit- 
tle (or no) injection of electromagnetically interacting en- 
ergy (i.e. high-energy 7's and e^'s). This may be seen 
from Fig.[5l where for larger life times Yx ^ 10~^^ is ruled 
out at the > 99% level. Here models are ruled out princi- 
pally due to violating the '^He/^H upper limit due to ''He 
photodisintegration. This fact is essentially not changed 
by the X-particle being charged, as may be seen by in- 
spection of the results in Ref. [14j for neutral decaying 
particles or by the dashed lines in Fig. [51 which shows 
the analogous limit for neutral particles. Immediately 
below the constraint line already between 20 — 80% of all 
models yield acceptable abundance yields. This param- 
eter space, is thus, currently not ruled out. However, it 
is conceivable that constraints on CHAMPs may be sig- 
nificantly tighted in future, to values possibly as low as 
Yx 'Z lO""'^ — lO""'*, in case rates for the important leac- 
tions'^EeCE-X- , X-fU, (^H-^-) + 2H-> (^H-X") 
+ iH, and CR-X-) + ^He^ (-^He-X") are deter- 
mined, and contrive to yield unacceptably large ^Li at 
low Yx ■ Here it is noted that the by far most important 
^Li producing reaction, potentially leading to such strin- 
gent constraints, is ^}Ie(^K—X~ , X~)^Li rather than 
2H(4He-X-,X-)6Li, yielding the bulk of the ^Li at 
temperatures T ^ 1 keV and not T « 6 — 8 keV. Impor- 
tant here is the CHAMP exchange reaction (^H— AT^) 
+ ^H^ (^H— A~) -|- ^H which may continously produce 
(^H— A~) bound states, (cf. Ref. [HI for more detail). 
This trend may be seen in Fig. [SI where ^Li/^Li isocon- 
tours are shown for a CHAMP with Bh = 0, and when 
the Born approximation for all rates is utilised. It is seen 
that for Tx ^ lO^sec factor ^ 1000 more ^Li is synthezised 
than for ^ lO^sec. 

It is evident from Figs. [4] and [5] that the probability 
distributions only flattens significantly for CHAMP life 
times Tx^b X lO^sec. For lO^sec ^Tx^ 5 x lO^sec the 
probability passes from > 99% observationally accept- 
able models to < 1% observationally acceptable models 
within one decade of Y^. This is because for shorter Tx 
possible late-time destruction of ^Li and ^Li may not yet 
be efficient due to large (^H— A^) bound state fractions 
only forming aXT"^! keV. Up to Tj, ^ 5 x lO^sec one may 
therefore approximately use the bounds on CHAMPS as 
given in Section 2. When Tx^5 x lO^sec late-time pro- 
cessing (destruction) of ^Li, ^Li, and ^H may occur ef- 



ficiently. For long life times it is thus proposed to use 
only the constraint imposed by possible ^He/^H overpro- 
duction. This constraint is not affected by bound states 
and the uncertainties in reactions including bound states. 
The ■^He isotope is special in that it may not be destroyed 
without destroying ^H as well. This is because the reac- 
tion •^He(^H— X~, A~)''Li is endothermic and other re- 
actions not involving ^H nuclei are Coulomb supressed. 
It is thus conservative to use the '^He/^H constraint, in 
particular, since there are other ^H destroying (but not 
producing) reactions within bound state nucleosynthe- 
sis. This translates into using the constraints given in 
Ref. for Bh = (and shown by the dashed curve for 
Tx ^ lO^sec in Fig. [5]). Since the fraction of rest mass Jem 
which is converted to electromagnetic interacting energy 
has been assumed /em = 1 for Bh = in the figures of 
Ref. [l3|, the constraint has to be rescaled accordingly, 
when neutrino losses are significant, or close-to mass de- 
generacy between mother and daughter particle exists. 
Finally, in the window 5 x lO^sec 'ZTx^i x lO^sec one 
may apply an additional constraint, stronger than that 
due to ■^He/^H overproduction. From Figs. [3| and [5| it 
is found that for > 3 x 10"" - lO'^^ only < 1% of all 
models are observationally acceptable. 

A procedure very similar to this has been very re- 
cently applied by Ref. to derive a lower limit on the 
gaugino mass parameter mi/2 in the constrained mini- 
mal supersymmetric standard model (CMSSM) when the 
gravitino is the lightest supersymmetric particle (LSP). 
This study considered only the catalyzed reaction ^H -|- 
(^He-A^) ^Li +X^ for intermediate life times of the 
next-to-LSP (NLSP), arguing that even without the ^Li 
bound long NLSP life times were already ruled out pri- 
orly by the ^He/^H upper limit. Since in the parameter 
space under investigation the NLSP is the stau, which 
has small hadronic branching ratio, the conclusion for 
Tx^5x lO^sec is thus not expected to change much when 
the present results are used. Using the results of the 
Monte-Carlo analysis shown in Figs. [4] and [5] typical stau- 
to-entropy ratios Yf ^ lO"'^'^ should be also disallowed for 
5 X lO^sec ^Tx^3x lO^sec due to ^Li overproduction and 
for Tj; ^ 3 X lO^sec due to '^He/^H overproduction, render- 
ing the conclusions of Ref. likely unchanged. 



IV. CONCLUSIONS 

In this letter the results of a detailed study of con- 
straints on charged massive particles X~ from Big Bang 
nucleosynthesis was presented. It was pointed out earlier 
that bound states between ^He and negatively charged 
X~ may lead to the efficient catalytic production of ^Li 
at T « 6 - 8keV Recently I have shown that 
BBN with CHAMPs enters a second phase of nucleosyn- 
thesis at T ^ 1 keV, capable of destroying all priorly syn- 
thesized ^Li. Altogether nineteen reactions important for 
late-time BBN were identified. When these processes are 
included in the analysis, drastic changes concerning lim- 
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its on the existence of CHAMPs when compared to those 
priorly derived 0, 0, 1 [M [H , 

are obtained. 

Limits on the existence of CHAMPs during and af- 
ter BBN are derived in two different decay time ranges. 
For very short Hfe times x lO^sec hmits are inde- 

pendent of the decaying particle being charged, or not. 
In the range 3 x lO^sec ^Tx^5 x lO^sec the important 
rates (i.e. ^HCHc— A"~, A"~)^Li) are relatively well de- 
termined such that a Monte-Carlo analysis may be 
avoided. It is stressed that the approximation of the 
^He-bound state fraction by the Saha equation leads to 
a factor 10 overestimate in the synthesized ^Li abun- 
dance. Due to substantial reaction rate uncertainties a 
full Monte-Carlo analysis had to be performed to obtain 
reliable and conservative bounds in the decay time range 
5 X lO^sec ^Tx^ lO^^sec. It was found that when a num- 
ber of reaction rates are large, and when electromagnetic 
energy injection is absent, baryon-to-cntropy ratios 
as large as ^ 10~^^ may be observationally acceptable. 
On the other hand, in case a number of reaction rates 
are determined more precisely, in particular the rates for 
4He(2H-X-,X-)6Li, (1H-X-) + (^H-X") + ^H, 
and (iH-X-) + ^He^ (iH-X") + ^He, it may be con- 
ceivable that limits on long-lived CHAMPs are improved 
by orders of magnitude, with the CHAMP-to-entropy ra- 
tio possibly constrained to be below 10""'^'^ — 10""'^*. A 
prescription is given for how to place conservative limits 
on CHAMPs, given current reaction rate uncertainties. 

The final result of this study is somewhat surpris- 
ing. When electromagnetic- and hadronic- energy release 
are included, and within the reaction rate uncertainties, 
conservative limits on charged decaying particles are no 
stronger than those on neutral particles. The exception 
here is the decay time range 10^— 3x lO^sec but only when 
the hadronic branching ratio is small Bh ^ 10~^ — 10~^, 
as is the case for supersymmetric staus. 

I acknowledge helpful discussions with and M. As- 
plund, S. Bailly, O. Kartavtsev, K. Kohri, A. Korn, 
G. Moultaka, M. Pospelov, J. Rafelski, G. Starkman, 
F. Steffen, V. Tatischeff, and T. Yanagida. 



V. ERRATUM: BOUNDS ON LONG-LIVED 
CHARGED MASSIVE PARTICLES FROM BIG 
BANG NUCLEOSYNTHESIS 



ring X~ from p to ^He, but also due to Coulomb barriers 
between p and nuclei being only partially shielded when 
protons are in bound states. The Monte Carlo analysis as 
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FIG. 7: Limits on the primordial CHAMP-to-entropy ratio 
Yx — nx/s (with nx-/s = Yx/2) as a function of their 
life time t^. Shown are constraint lines for CHAMPs of 
mass Mx = 1 TeV and a variety of hadronic branching ra- 
tios Bh ~ 10"^ — 1, as labeled in the figure. Solid (red) lines 
correspond to the conservative limit ®Li/^Li < 0.66, whereas 
dashed (blue) lines correspond to ^Li/'^Li < 0.1. It is seen 
that only for CHAMPs with Bh < 10"^ the effects of bound 
states become important. For small decay times Tx ^ 10'^ sec 
and large decay times Tx ^ lO^sec the limits on CHAMP abun- 
dances are virtually identical to those on the abundance of 
neutral relic decaying particles [l^. 



X 




A recent detailed computation by Kamimura et al. [20j 
of the most important reaction rates entering bound state 
nucleosynthesis (see also the discussion in ReL [2l|) has 
established that late time r ^ lO^sec processing of light 
element abundances as envisioned possible by the present 
author [l^l does not usually take place. This is mostly 
because of the abundance of proton- A ~ bound states 
staying small due to efficient exchange reactions transfer- 



FIG. 8: As Fig. [3 but for Mx = 100 GeV. 



performed in Section II and III is therefore superfluous. 
Figs. 4-6 of the paper are replaced by Figs. 7 and 8 in 
this erratum. It is cautioned, however, that for relatively 
large Yx > 2l4jjc late time processing may still occur to 
some degree. 
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